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ABSTRACT 



The surface gravities of cool (T^ff < 13, 000 K) hydrogen-atmosphere DA white dwarfs, determined from spectroscopic analyses, are 
found to be significantly higher than the canonical value of log ~ 8 expected for these stars. It was recently concluded that a problem 
with the treatment of convective energy transport within the framework of the mixing-length theory was the most plausible explanation 
for this high-log g problem. We pursue the investigation of this discrepancy by computing model spectra of cool convective white 
dwarfs from a small sequence (1 1,300 K < T^ff < 12,800 K) of 3D hydrodynamical model atmospheres, which feature a sophisticated 
treatment of convection and radiative transfer. Our approach is to proceed with a differential analysis between 3D and standard ID 
models. We find that the 3D spectra predict significantly lower surface gravities, with corrections of the right amplitude as a function 
of efi"ective temperature to obtain values of log g ~ 8 on average. We conclude that the surface gravity distribution of cool convective 
DA white dwarfs is much closer to that of hotter radiative objects when using, for the treatment of the convection, 3D models instead 
of the mixing-length framework. 
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1. Introduction 

The shape of Balmer lines in hydrogen-atmosphere DA white 
dwarfs is particularly sensitive to variations of the atmospheric 
parameters (Tgff and \ogg). Therefore, the spectroscopic tech- 
nique, which consists in comparing the observed line profiles 
of the Balmer series with the predictions of detailed model at- 
mospheres, is by far the most accurate method for determining 
the atmospheric parameters. This technique was initially applied 
to a sample of 37 coo l (Teff < 13,000 K) DA white dwarfs by 
[Bergeron et al.l ([1990') who showed that the log g values mea- 
sured for these white dwarfs were significantly higher than the 
expected canonical value of log g ~ ^, which is about the mean 
value determined later for hotter DA stars (e.g., [Bergeron et al.l 
Il992h . This discrepancy is now observe d in all large spectro- 
scopic surveys of DA white dwarfs (e.g., Eisenstein et al.ll2b06t 
iGianninas et al.ll2"009h and no satisfactory explanation has been 
reported until now. In Fig. [1] we present our current view of 
this high-log 9 problem for the spectroscopically identified DA 
stars in the Sloan Digital Sky Survey (SDSS). This long-standing 
problem implies that we cannot obtain reliable atmospheric pa- 
rameters determined from spectroscopy at the cool end of the 
white dwarf evolutionary sequence. Since the problem manifests 
itself for white dwarfs older than ~1 Gyr, it has important impli- 
cations on our ability to use these objects as cosmochronome- 
ters a nd distance indicators dWinget et al.lll987l : [Fontaine et al.l 
I2001h . 




Te„ (K) 

Fig. 1. Surface gr avity distribution as a function of T^ff for the 
SDSS sample; see'Tremblav et al.l (1201 Ih for details on these de- 
terminations. An evolutionary model fromfFontaine et al. (200i]) 
at the median mass of the sample (0.59 Mq) is shown as a dashed 
line. 



Since white dwarf stars are expected to cool at constant 
mass and almost constant radius, the sudden and significant in- 
crease in logg can not be explained in terms of simple astro- 
phy sical argument s. Tremblav et alj (j2010, hereafter TBIO; see 
also iKoester et al.ll2009h made an extensive review of the possi- 
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ble solutions, which range from inadequate assumptions about 
the composition of these stars to inaccuracies in the model at- 
mosphere calculations. The non-detection of He i lines in high- 
resolution Keck observations of cool DA stars by TBIO ruled 
out the systematic presence of helium in the atmospheres, which 
would mimic high er spectroscopi c log g determinations, as ini- 
tially proposed by Ber geron et alJ (IT990). TBIO concluded - as 
[Koester et al. (2009) did - that a problem with the treatment 
of convective energ y transport, currently based on the mixing- 
length theory (MLT: lB5hm-Vitens3l 1 958h . is the only viable ex- 
planation for the high-log g problem. 

The MLT approximation is a phenomenological description 
of the complex convective fluid movements in a star According 
to this model, bubbles transport energy in the convection zone 
over a characteristic distance - the mixing length - where they 
dissolve and release their extra energy in the environment. Even 
though the MLT framework, properly calibrated, provides a rea- 
sonable description of white dwarf model atmospheres (TBIO 
and references therein), even small changes in the convective ef- 
ficiency, or equivalently the mixing-length parameter, can have a 
significant impact on the predicted line profiles (see, e.g.. Fig. 3 
of lBergeron et alj 19951) in the regime where the high-log g prob- 
lem is most important. This comes from the fact that convec- 
tive energy transport may affect dramatically the temperature 
and pressure structures of the atmosphere, and consequently the 
shape of the broad hydrogen lines, formed at different depths 
in the atmosphere. The observations reveal that the high-log 
problem seems to be related to the line profiles only, and not to 
the continuum flux (as demonstrat ed from photometric analyses 
of these stars; see, e.g., TBIO and lKoesteret"ani200 9^. and that 
the problem occurs when convection becomes important in the 
photosphere (TBIO). This is consistent with the MLT framework 
being the source of the problem. 

A more detailed multidimensional radiation-hydrodynamics 
(RHD) treatment of convective motions in DA white dwarfs, 
ba sed on first physica l principles, h as already been studied 
bvlLudwig et alj (ll994 . ISteflren et al.1 d 19951) . and lFrevtag eTal] 
(Il996h . These studies revealed that the diff'erences between 2D 
spectra and those computed from standard ID hydrostatic mod- 
els were qualitatively small. However, a visual inspection of 
Fig. 13 from TBIO reveals that the model predictions do not need 
to change by much to rectify the high-log g problem. These 2D 
calculations have not been pursued further, and in the meantime, 
full three dimensional RHD simulations (hereafter 3D models) 
have become more sophisticate d, and accurate enou gh to make 
realistic predictions (see, e.g., lAsplund et alj 120001 and refer- 
ences therein). 

In this work, we present an updated investigation of the dif- 
ferences between standard DA spectra and those calculated from 
3D model atmospheres. In Sect. 2, we describe our computation 
of a small number of 3D model atmospheres for cool convec- 
tive DA white dwarfs. We then address, in Sect. 3, the status 
of the high-log g problem from a purely differential analysis be- 
tween the Balmer line profiles predicted by the 3D models, and 
those computed from ID models with the same underlying mi- 
crophysics. 



2. Improved 3D model atmospheres 

We computed a smaU set (11,300 K < T^s < 12,800 K) of 
pure-hydrogen 3D model atmospheres at log g - % with th e 
CO^BOLD codfl (iFrevtag et alJl2002l: IWedemever et alJl2064l) . 



This code solves the time-dependent hydrodynamical equations 
(conservation of mass, momentum, and energy) for a fully com- 
pressible fluid with a finite-volume approach, coupled to the ra- 
diative transfer The DA white dwarfs investigated here have 
a shallow convective zone, which can be included entirely in 
the simulation domain (for the vertical direction). This allows 
for dynamically stable layers at the top and bottom of the do- 
main, and accounts for the entire overshooting layer The lateral 
boundaries are periodic, and the radiative flux entering the do- 
main at the bottom defines the effective temperature. We rely 
on a standard stellar equation of state, and the opacities are th e 
same as those used in the simulations of Ludwig et al.l (|1994|) . 
We computed both gray and non-gray m odels, the latter using 
7 opacity bins for the radiative transfer (iLudwig et alj|1994l) . 
Table [1] presents the basic properties for our set of models. 
In Fig. [2] we illustrate the structure of the non-gray 11,975 K 
model. We find th at our new models ar e qualitatively similar to 
those computed bv lLudwig et all (11994) and Steffen et al.. (.19951) 
in terms of the properties of the convective zone (i.e., depth and 
maximum velocities). 

Table 1. Parameters of our 3D model atmospheres 





logf 


Opacity bins 


Time" 




(K) 






(s) 


(%) 


Non-gray models 


11,300 


8.0 


7 


10 


2.46 


11,975 


8.0 


7 


10 


0.07 


12,390 


8.0 


7 


10 


0.16 


12,800 


8.0 


7 


10 


0.02 


Gray models 


11,350 


8.0 


1 


40 


1.20 


11,975 


8.0 


1 


10 


0.70 


12,390 


8.0 


1 


15 


0.03 


12,800 


8.0 


1 


15 


0.06 



www.astro.uu.se/~bf/co5bold_main.html 



Total stellar time. Maximum flux error, defined as (max|F(z)/F* - 
1|), taken over the whole depth of the computation box. F' is the nom- 
inal flux crT^g, and F{z) is the total energy flux (radiative, convective, 
kinetic, viscous and potential fluxes) averaged over 2 seconds (or 5 sec- 
onds for the coolest gray and non-gray models) at the end of the simu- 
lation. Note that the maximum flux error generally occurs at the base of 
the simulation and that in the line forming layers, the flux error is much 
smaller. 

Notes. The box size is 7.45 km (150 grid points) for both horizontal 
directions and 4. 17 km for the vertical direction (100 grid points). 



In order to make our comparison strictly differential, we 
also employ ID models in which the microphysics and radiation 
transfer numerical schemes (e.g., the opacity binning) are iden- 
tical to those included in the 3D models. W e computed these 
hydro static ID models with the LHD code ( ICaffau & Ludwigl 
l2007h . which treats convection within the MLT framework us- 
ing our most recent ML2/a - 0.8 parameterization (see TBIO 
for details). In Fig.[3j we present these models hereafter named 
IDmlt (the gray models are omitted for clarity). In addition, 
we computed mean ID structures that correspond to the (r4)i/4 
and (P) spatial average of the 3D models over surfaces of con- 
stant Rosseland optical depth. We followed this procedure with 
temporal averages over 12 selected snapshots. These averages, 
called (3D) models and shown in Fig.[3j have had their horizon- 
tal temperature and pressure fluctuations removed, and therefore. 
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Fig. 3. Temperatures versus Rosseland optical depth for (3D) 
(solid lines) and IDmlt (dashed lines) non-gray simulations at 
Teff = 11,300, 1 1,975, 12,390 and 12,800 K (shifted on the ver- 
tical axis by 0, 2, 4 and 6 kK, respectively, for clarity). 



a comparison of 3D and (3D) models highlights the effect of 
these fluctuations. 

At first, we rely on the Linfor3D three-dimensional spec- 
tral synthesis codq3 to compute H/3 line profiles from our DA 
model atmospheres. Figure |4]presents the 3D, (3D), and IDmlt 
profiles at T^s = 11,300 K and 12,390 K. It is already obvi- 
ous that the 3D spectra are indistinguishable from the spectra 
computed from the corresponding averaged (3D) models, an ef- 
fect also observed by Steffen et al. (1995) in their 2D simula- 
tions. This is a non-trivial result since significant temperature 
and pressure inhomogeneities are present in the hydrodynamical 
simulations (see Fig. |2]l. In the case of main-sequence stars, for 
instance, differences between 3D and (3D) spectra can be im- 
portant (Ludwig et al. 2009). However, in the case of DA white 
dwarfs, an optical spectrum can be safely computed using a 
single averaged atmospheric structure. In this context, one can 
therefore neglect the time-consuming task of 3D spectral syn- 
thesis. 

On the other hand. Fig. |4] also demonstrates that the differ- 
ences between the (3D) and IDmlt models are significant. This 
suggests that most of the effects from the RHD treatment of con- 
vection originate from changes in the mean vertical atmospheric 
structure. An examination of Fig.|3]indeed reveals that the (3D) - 
hereafter equivalent to 3D - and IDmlt structures have a signif- 
icantly different shape in the photosphere. For the remainder of 
this study, our goal is to quantify the differences between these 
models. 

We chose to compute optical synthetic spectra (3500 A < 
X < 5500 A) for all our 3D and IDmlt structures calculated 
in this work using the standard spectral sy nthesis code for DA 
white dwarfs from the Montreal group (Tremblay & Bergeron 
l2009i. TBIO). The rationale behind this choice is that, as de- 

2 www.aip.de/~mst/linfor3D_main.html 
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Fig. 4. Comparison of H/3 line profiles calculated with Linfor3D 
for the 3D (solid lines, averaged over 12 selected snapshots) and 
IDmlt (dashed lines) models in the case of two non-gray sim- 
ulations at Teff of 11,300 K and 12,390 K (the latter shifted by 
0.2 units of flux for clarity). The dotted lines (indistinguishable 
from the 3D spectra in the figure) also show the (3D) profiles re- 
sulting from spatial and temporal averages of the 3D structures. 



scribed above, one can neglect the 3D spectral synthesis, and 
moreover, our code relies on an equation of state and opacity 
sources that are more appropriate for DA white dwarfs (e.g., 
with a proper account of line blending, including non-ideal gas 
effects). We verified that the structures calculated in this work 
from the LHD and CO^BOLD codes (with the microphysics dat- 
ing from lLudwig et al.]|I994l) are comparable to those computed 
with the Montreal code. This implies that the differences in the 
microphysics are small, after all, and that our procedure is ade- 
quate for a differential analysis. 

3. Application to the high-log ^ problem 

To quantify the effects of the 3D white dwarf models, we com- 
pute in this section the 3D atmospheric parameter corrections 
with respect to the reference IDmlt models. These are defined as 
the differences in Teff and log g when we fit the normalized line 
profiles of both the IDmlt and 3D spectra with a reference gricQ. 
More specifically, the corrections are defined as 1Dmlt-3D, 
since we are interested in how much the atmospheric parameters 
of real stars would change w hen using 3D model spectra. Our 
approach is similar to that of iLudwig et alj (l2009l) who studied 
these coiTections for late-type main sequence stars. We use as 
a reference grid the spectra computed in TBICQ. Obviously, the 
atmospheric parameter corrections will not be exact since the 

^ We neglect the center of the lines (\AA\ < 3 A), which are formed at 
low values of Tr where the number of frequencies taken into account in 
the structure calculations can have a strong impact on the predictions. 

* We also calculated a similar reference grid from gray structures to 
fit our gray models. 
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independent reference grid does not exactly reproduce the spec- 
tra computed in this analysis. Nevertheless, the good quality of 
our fits achieved in all cases implies that our approach should 
be valid to within a few percents. The 3D log g corrections are 
presented in Fig.|5] In the case of the 3D corrections in effective 
temperature, we find that they are all relatively small (~ 200 K 
at Teff ~ 12, 000 K). They are certainly of the same order as the 
accuracy of the spectroscopic method, and therefore we do not 
discuss them further We note that the 3D corrections should ap- 
proach zero for T^s > 13,000 K, since the convection zone is 
rapidly shrinking, and the IDmlt and 3D spectra are otherwise 
identical. 



0.00 I 




13000 12500 12000 11500 11000 

Te„ (K) 

Fig. 5. Non-gray (filled circles) and gray (filled squares) 3D log g 
corrections. In comparison, we show the shifts in surface gravity 
required to obtain a stable log g distribution as a function of Teff 
using the samples of the SDSS (open circles, also displayed in 
Fig.lTJ and Gianninas et al. (2009; open squares). The points are 
connected for clarity. 



The non-gray 3D log g corrections are the key values to as- 
sess the high-log g problem. A preliminary observation reveals 
that these corrections are negative for all spectra (i.e., the 3D 
models predict lower surface gravities), which is in the appro- 
priate direction to correct the problem. The mean amplitude of 
the corrections is significant (~ -0.15 dex), which may not be 
surprising since the complicated flow patterns are not expected 
to be fully represented by a unique free parameter for all atmo- 
spheric depths, and for all white dwarfs. Also given in Fig. |5] 
are the observed shifts in \ogg derived from the SDSS white 
dwarf sample displayed in Fig.lT] These shifts correspond to the 
corrections required, in a bin of 1000 K around the model tem- 
perature, to match the median mass value obtained from hot DA 
stars, represented by the dashed line i n Fig. [1] We also com - 
puted similar shifts from the sample of iGianninas et al.l (l2009h . 
drawn from the Villanova White Dwarf Catalog, to account for 
observational uncertainties. We can see that the 3D corrections 
have almost exactly the right slope and amplitude to solve the 
high-log g problem. This result clearly indicates that the surface 
gravity distribution of cool convective DA stars (Teff > 11,000 



K) would be much more stable, as a function of Tgff , if 3D model 
atmospheres were used instead of the ID standard models. 

We also verified that the non-gray log g corrections are fairly 
independent of the Teff corrections by redoing our fits at fixed 
temperatures (Teff - T^d)- The corrections we obtain with this 
method are within ~ 20% of those found in Fig. |5] Finally, the 
fact that we have a set of 3D gray models allows us to evaluate 
the impact on our results of the opacity binning procedure in the 
non-gray models. The shifts in log g that we obtain with 1 opac- 
ity bin (gray models) are smaller than those with 7 bins (non- 
gray models), but still correct to a large degree the high-log gr 
problem. This shows that further improvements in the binning 
procedure of the RHD models are not expected to change our 
differential analysis significantly. 

4. Conclusion 

We conducted a differential analysis of synthetic spectra for cool 
convective DA white dwarfs computed from model atmospheres 
with a detailed 3D radiation-hydrodynamical treatment of con- 
vective motions, and with the standard ID mixing-length theory. 
We introduced the concept of 3D atmospheric parameter correc- 
tions, which are defined as the difference in Teff and log g be- 
tween IDmlt and 3D spectra when fitted with an independent 
grid. These corrections illustrate by how much the atmospheric 
parameters of DA stars would change if we were to use a grid 
of 3D spectra to fit the observations. We find that the 3D surface 
gravity corrections have the correct amplitude (~ 0.15 dex) to 
solve the high-log g problem for DA white dwarfs in the range 
11,300 K < Teff < 12,800 K. This indubitably suggests that a 
weakness with the standard MLT theory is creating the high log g 
values reported in recent spectroscopic analyses (e.g., TBIO). 

Our next goal will be to review the microphysics and the 
opacity binning procedure in the CO^BOLD code to connect the 
absolute properties of these hydrodynamical models with those 
of standard white dwarf models. We will then be able to compute 
improved 3D model spectra that can actually be compared with 
observations to determine the atmospheric parameters. We also 
plan to look at convective white dwarfs cooler than those studied 
in this work, where the MLT theory predicts an increasingly adi- 
abatic convection, which is less sensitive to the mixing-length 
parameter. It remains to be seen how the 3D log g corrections 
would behave in this regime. 
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X [km] 

Fig. 2. Temperature structure for the 3D white dwarf model at Teff = 11,975 K and \ogg - 8, for a slice in the horizontal-vertical 
xz plane through a box with coordinates jc, t/, z (in km). The temperature is color coded from 60,000 K (yellow) to 8000 K (blue). 
The arrows represent relative convective velocities (integration time of 20 ms), while thick lines correspond to contours of constant 
Rosseland optical depth. 



